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Abstract. We present preliminary numerical evidence that the physical con- 
ditions in high-mass star forming regions can arise from global gravitational 
infall, with the velocity dispersions being caused primarily by infall motions 
rather than random turbulence. To this end, we study the clumps and cores 
appearing in the region of central collapse in a numerical simulation of the for- 
mation, evolution, and subsequent collapse of a dense cloud out of a transonic 
compression in the diffuse atomic ISM. The clumps have sizes ~ 1 pc, masses of 
several hundred Mq, and three-dimensional velocity dispersions ~ 3 km s" 1 , in 
agreement with typical observed values for such structures. The clumps break 
down into massive cores of sizes ~ 0.1 pc, densities ~ 10 5 , masses 2-300 M©, 
with distributions of these quantities that peak at the same values as the mas- 
sive core sample in a recent survey of the Cygnus X molecular cloud complex. 
Although preliminary, these results suggest that high-mass star forming clumps 
may be in a state of global gravitational collapse rather than in equilibrium 
supported by strong turbulence. 



1. Introduction 

High-mass star forming regions are character ized by more extreme physical con- 



riign-mass star tormmg regions are cnaracter izea by m ore extreme pnysica l con- 
ditions than their low-mas s counterparts (e.g. iGaray Sz Lizanol 11993 : iKurtz et al 



200d : iBeuther et alJl2007T ). having "clumps" of sizes 0.2-0.5 pc, mean densities 



n ~ 10 cm , masses between 100 and 1000 Mq, and velocity dispersions rang- 
ing between 1.5 and 4 km s . In turn, these clumps break down into even denser 
"cores" that are believed to be the immediate precursors of single or gravitation- 
ally bound multiple massive protostars. The high velocity dispersions of these 
clumps are generally interpreted as strong turbulence that manages to suppor t 
the clumps against gravity (e.g., iGarav Lizanol Il999l : iMcKee Tanl feOQS). 
However, the notion of "turbulent support" is difficult to maintain at the scales 
of these cores. Turbulence is a flow regime in whic h the largest velocity differ- 
ences are associated with the largest separations (e.g. lFriscn1ll995l ). and moreover 
in the case of supersonic turbulence the clumps are expected to be formed by 
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large-scale compressive motions, so that the turbulence is likely to have a strong 
compressive component dHunter k Flecklll982l ; iBallesteros-Paredes et al.lll999al . 
l2008l : IVazauez-Semadeni et al.ll2008l ~ So, it is difficult to imagine a turbulent ve- 
locity field inside the clumps that is completely random in such a way as to only 
provide support against gravity - the compressive component may rather foster 
core contraction. In ad dition, numerical simulations of cloud formation in the 



1995, 



19961: iPassot et al.lll995l: 



diffuse atomic ISM (e.g. ly azque z-Semadeni et al .._ 
Ballesteros-Paredes et al.lll999alJbl: lAudit Hennebelle 1120051: iHeitsch et al.ll2005l . 
20061 : IVazauez-Semadeni et alJl2006l . l2007l : lHennebelle fc AuditfeOOTf) and of star 



form ation in turbulent, self-gravitating c l ouds fe.g.lKlessen. Heitsch &: Mac Low 

20001 : IHeitsch. Mac Low fe Klessenll200ll : iBate. Bonnell. fc Brommll2003l : IVazquez-Semadeni et al 



20051.12007) show that the velocity fields are organized at all scales, exhibiting a 



continuity from the large scales outside the clumps all the way to their interiors. 

The large-scale compressions can be of turbulent origin (e.g., passing spi- 
ral arm shocks, supernova shells, or simply the general transonic turbulence 
in the diffuse medium) or of gravitati onal origin (e.g., l arge-scale gravitational 
or m a gneto-gravitational instab i lities: lElmegreen 1991 : Kim. Ostriker Sz Stone 
20021 : Field. Blackman. Ketol I2008T ) . The self-gravitating simulations men- 



tioned abo ve exhibit gravitationally d riven motions up to the largest scales. In 
particular, IVazquez-Semadeni et all ( 20071 . hereafter Paper I) have presented 
simulations of molecular cloud formation by generic compressions in the dif- 
fuse atomic ISM and of its sub sequent collapse and star-forming stage, using the 
SPH/N-body code GADGET (jSpringel et a,l.ll20Qlh . complemented with param- 
eterized heating and cooling taken fromlKovama Inutsuka ( 20021 ) and a sink- 
particle prescription (jjappsen et al.ll2005l ). In these simulations, the clouds that 
formed acquired their initial turbulence from instabilities of the compressed layer 
(IVishniadll99llWalder fc Folinilll998|:lKoYama & Inutsuka 1l2002l: IHei tsch et al 
120051 . 120061 : IVazquez-Semadeni et ali i 
ally unstable and began contracting. 



but soon they became gravitation- 
This contraction phase was nevertheless 



characterized by a virial-like energy balance with |-E gr av| ~ ^E^ D , which was 
however due to the gravitational contraction, not to virial equilibrium. During 
the global contraction, clumps produced by the initial turbulence proceeded to 
collapse on their own, forming what resembled low-mass star forming regions. 

In Paper I, we speculated that the global contraction might be halted by 
stellar energy feedback before the global collapse was completed. However, here 
we forgo that speculation, and take the simulation at face value, presenting 
a preliminary study of the physical conditions in the region where the global 
collapse finally converges, showing that they resemble the physical conditions 
of high-mass star forming regions, thus suggesting that such regions may be in 
generalized gravitational collapse rather than in a state of turbulent "support" . 



2. The numerical model 

The simulation we consider is the one labeled L256Af0.17 in Paper I. We refer 
the reader to that paper for details. Here we just mention that it is an SPH simu- 
lation with self-gravity, parameterized heating and cooling implying a thermally 
bistable medium, using 3.24 x 10 6 particles, using sink particles, and initially 
set up to produce a collision of streams of diffuse gas (at the same density as 
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Figure 1. Column density plot of the central 50 pc of this simulation in 
the y-z plane at t = 23.4 Myr, integrating over the central 8 pc along the x 
direction. The box shows the region analyzed in |J3T] The dots show the stellar 
objects (sink particles) already formed in the simulation by this time. Radial 
filamentary streams of intermediate-density gas are seen to still be accreting 
onto the central cloud. 



their surroundings) that induces a transition to the cold, dense phase and the 
generation of turbulence in the dense gas. The size of the numerical box was 
256 pc, and the inflow velocity of the colliding streams was 1.25 times the sound 
speed in the ambient gas, which had T = 5000 K and a mean density of n = 1 
cm" 3 . The turbulent cloud eventually reached densities typical of molecular gas, 
and began to contract gravitationally. The cloud had a flattened shape during 
most of its evolution, and star formation began at t ~ 17 Myr in the periphery 
of the cloud, where secondary compression produced by the gas squirting off 
the collision site produced the highest initial densities. Animations showing the 
large-scale evolution of this simulation can be found in the electronic edition of 
Paper I. 

By t = 23.4 Myr, the global collapse is completed, although the residual 
turbulence causes the motions to have a random component, so that the collapse 
center spans several parsecs across. Figure [T] shows a column density of the cen- 
tral 50 pc of this simulation in the y-z plane at t = 23.4 Myr, integrating over 
the central 8 pc along the x direction. The dense cloud is seen near the center of 
the image, with streams of gas still infalling onto it. A region 8 pc on a side con- 
taining the cloud is indicated by the square, for which animations can be found at 
: //www . astrosmo . unam . mx/ . vazquez/turbulence/movies . html[ These 



Figure 2. Two views of the central 8-pc cubic region highlighted in Fig. 
[U Left panel: Column density integrated along the z direction. Right panel: 
Column density integrated along the y direction. The two squares show the 
regions called Clump A and Clump B. 



animations show the evolution of the central 8 pc of the simulation for 22.1 < 
t < 24.7 Myr, showing a violent collapsing evolution in which inf ailing clumps 
of gas interact but do not entirely merge, but rather undergo shredding and 
distortion. It is this violent region that we analyze in the next section. 



3. Physical conditions of clumps and cores in the collapse center 

Figures [2^ and [2)3 show two views of two parsec-sized clumps within the 8-pc 
cloud, to we refer as "Clump A" and "Clump B". Note that these are just cubic 
boxes enclosing the dense clumps, rather than actual clumps defined by any 
clump-finding algorithm. Their properties can be compared with the observed 
typical properties of clumps in high-mass star forming regions. We do this by 
interpolating the SPH data for the central 8-pc region into a fixed grid with a 
resolution of 256 3 . Clump A, which has a size of 1.5 pc per side, has a mass 
M = 1400Mq, a mean density of (n) = 1.27 x 10 4 cm -3 , and a three-dimensional 
velocity dispersion a = 3.6 km s _1 . Clump B, in turn, has a linear size of 0.8 
pc, a mass M = 300MQ, a mean density of (n) = 1.72 x 10 4 cm -3 , and a 
velocity dispersion a = 2.8 km s~ 4 . So, in general these properties compare well 
with those quoted in §Lj except perhaps for slightly lower mean densities than 
typical, which can be understood as a consequence of our usage of cubic boxes 
rather than clumps. The boxes include some lower-density gas. On the other 
hand, our densities fare in well with th ose reported for ma ssive starless clumps 
in the Cygnus X molecular complex bv lMotte et al.l ( 20071 ). 



Within these clumps, we identify cores by applying a simple clump-finding 
algorithm based on finding connected sets of grid points whose densities are 
above a certain threshold. For this preliminary analysis, we consider a single 
threshold n t h r = 5 x 10 4 cm -3 , leaving us with 20 cores, 14 of which have 
M > 4Mq. For each core, we measure its mass, mean density, and velocity 
dispersion, and estimate its size as ~ (3V/47T) 1 / 3 , where V is its volume. 
These properties can be compared with those reported b y recent sur v eys o f 
cores in high-mass star forming regions, such as that by Motte et al. (2007) 
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Figure 3. Histograms of the size (top left), mass (top right), mean density 
(bottom left) and three-dimensional velocity dispersion (bottom right) of the 
cores within Clumps A and B of the simulation (solid lines), normalized to 
the peak of the histogram. Dotted lines: same for the core sample in Cygnus 
X reported bv iMotte et al.l (|2007t ). except for the velocity dispersion, since 
their observations were performed in the 1.2 mm continuum. The peaks of 
the distributions are seen to coincide for size, mass and mean density. The 
reasons for the simulation distributions being narrower are discussed in the 
text. 



for the Cygnus X region. These authors give these data for a sample of 129 
dense cores with sizes ~ 0.1 pc, masses 4-950 M@, and mean densities ~ 10 5 
cm -3 . In order to compare their data to the cores in our clumps, we fortuitously 
select the first 57 cores (those in the first page) in their Table 1 and show the 
distribution of their properties in the histograms presented in Fig.[3]by the dotted 
lines. Superimposed on these histograms, the solid lines show the corresponding 
distributions for the 14 cores more massive than 4Mq in Clumps A and B of 
our simulation. We see that, although the distributions of the Cygnus X cores 
are in general broader than those of the cores in our simulation, the peaks of 
the distributions match for size, mass, and mean density. 

The fact that the distribution of densities for the cores in the simulation 
is narrower than the Cygnus X one is most likely a result of our sample being 
significantly smaller than that for Cygnus X, and of our having considered only a 
single density threshold for defining the cores. Our previous experience with this 
clump-finding method is that the mean density of the cores found is generally 
within less than one order of magnitude of the threshold density. This limitation 
is also the probable cause of the absence of the small-size tail of the cores in the 
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simulation, since smaller cores should appear when higher density thresholds are 
considered. So, in order to obtain a wider range of densities, masses and sizes, 
it is necessary to consider a suite of density thresholds. We plan to do it in the 
final study, to be presented elsewhere. 



4. Conclusions 

In this contribution we have presented preliminary numerical evidence that the 
physical conditions in high-mass star forming regions can arise from global grav- 
itational infall, with the velocity dispersions being caused primarily by infall 
motions rather than random turbulence. The evidence comes from the first 
study of core properties in a simulation of the entire evolution of a molecular 
cloud, from its formation in the diffuse atomic ISM to its gravitational col- 
lapse. Although the analysis presented here is only preliminary, it is consistent 
with recent suggestions, based on comparisons between simula tions and observa- 
tions, that molecular clouds (e.g., lHartmann Sz BurkertJ 12007) and clumps (e.g. 
Peretto. Hennebelle Andre 120071 ) may be in a state of gravitational collapse. 



If confir med, these suggesti o ns po int towards a return to the original sugges- 
tion by iGoldreich fe Kwaii ( 19741 ) that the observed linewidths in molecular 



clouds are due primarily to grav i tation al contraction. This suggestion was dis- 
missed by IZuckerman fc Palmerl (|l974l ) through the argument that this would 
imply a much larger average star formation rate in the Galaxy than observed. 
However, it is possible that this criticism may be overcome if magnetic field 
fluctuations in the clouds imply that some parts of them are magnetically sup- 
ported so that only the non-supported parts parts of the clouds un dergo collapse 



(jHartmann. Ballesteros-Paredes. Sz Bergin Il200ll ; lElmegreenll2007l ) , and the star 



formation rate is then regulated by stellar feedback, with globally collapsing mo- 
tions arising only for those regions that manage to "percolate" through the field 
fluctuations and the stellar-feedback motions. We plan to perform simulations 
including magnetic support and stellar feedback in future studies to investigate 
this possibility. 
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